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Abstract 

Associated Legendre functions arise in many problems of mathematical physics. By using 
the generalized Abel-Plana formula, in this paper we derive a summation formula for the 
series over the zeros of the associated Legendre function of the first kind with respect to 
the degree. The summation formula for the series over the zeros of the Bessel function, 
previously discussed in the literature, is obtained as a limiting case. The Wightman function 
for a scalar field with general curvature coupling parameter is considered inside a spherical 
boundary on background of constant negative curvature space. The corresponding mode 
sum contains series over the zeros of the associated Legendre function. The application of 
the summation formula allows us to present the Wightman function in the form of the sum 
of two integrals. The first one corresponds to the Wightman function for the bulk geometry 
without boundaries and the second one is induced by the presence of the spherical shell. For 
points away from the boundary the latter is finite in the coincidence limit. In this way the 
renormalization of the vacuum expectation value of the field squared is reduced to that for 
the boundary- free part. 



PACS numbers: 02.30.Gp, 03.70.+k, 04.62,+v 



1 Introduction 

In a number of problems in mathematical physics we need to sum over the values of a cer- 
tain function at integer points, and then subtract the corresponding integral. In particular, in 
quantum field theory the expectation values for physical observables induced by the presence of 
boundaries are presented in the form of this difference. The corresponding sum and integral, 
taken separately, diverge and some physically motivated procedure to handle the finite result, 
is needed. For a number of boundary geometries one of the most convenient methods to obtain 
such renormalized values is based on the use of the Abel-Plana summation formula [TJ [2] (for 
different forms of this formula discussed in the literature see also [3]). Applications of the Abel- 
Plana formula in physical problems related to the Casimir effect for flat boundary geometries 
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and topologically non-trivial spaces with corresponding references can be found in [HE]. The 
use of this formula allows to extract in a cutoff independent way the Minkowski vacuum part 
and to obtain for the renormalized part rapidly convergent integrals useful, in particular, for 
numerical calculations. 

However, the applications of the Abel-Plana formula in its standard form are restricted 
to the problems where the normal modes are explicitly known. In [6] we have considered a 
generalization of this formula, which essentially enlarges the application range and allows to 
include problems where the eigenmodes are given implicitly as zeros of a given function. Well 
known examples of this kind are the boundary- value problems with spherical and cylindrical 
boundaries. The generalized Abel-Plana formula contains two meromorphic functions and by 
specifying one of them the Abel-Plana formula is obtained (for other generalizations of the 
Abel-Plana formula see [5j[71[8]). Applying the generalized formula to Bessel functions, in [HE] 
summation formulae are obtained for the series over the zeros of various combinations of these 
functions (for a review with physical applications see also Ref. [3j QUI E] ) • 

The summation formulae derived from the generalized Abel-Plana formula have been ap- 
plied for the evaluation of the vacuum expectation values of local physical observables in the 
Casimir effect (for the Casimir effect see [HO [12]) for plane boundaries with Robin or non-local 
boundary conditions [13J, for spherical boundaries in Minkowski and global monopole bulks [14] 
and for cylindrical boundaries in Minkowski and cosmic string bulks [15]. By making use of 
the generalized Abel-Plana formula, the vacuum expectation values of the field squared and 
the energy-momentum tensor in closely related but more complicated geometry of a wedge with 
cylindrical boundary are investigated in [16] for both scalar and electromagnetic fields. As in the 
case of the Abel-Plana formula, the use of the generalized formula in these problems allows to 
extract the contribution of the unbounded space and to present the boundary-induced parts in 
terms of exponentially converging integrals. In [T7] summation formulae for the series over the 
zeroes of the modified Bessel functions with an imaginary order are derived by using the gener- 
alized Abel-Plana formula. This type of series arise in the evaluation of the vacuum expectation 
values induced by plane boundaries uniformly accelerated through the Fulling-Rindler vacuum. 
Another class of problems where the application of the generalized Abel-Plana formula provides 
an efficient way for the evaluation of the vacuum expectation values is considered in [18] . In 
these papers braneworld models with two parallel branes on anti-de Sitter bulk are discussed. 
The corresponding mode-sums for physical observables bilinear in the field contain series over 
the zeroes of cylinder functions which are summarized by using the generalized Abel-Plana for- 
mula. The geometry of spherical branes in Rindler-like spacetimes is considered in |19] . In [20] 
from the generalized Abel-Plana formula a summation formula is derived over the eigenmodes 
of a dielectric cylinder and this formula is applied for the evaluation of the radiation intensity 
from a point charge orbiting along a helical trajectory inside the cylinder. 

The physical importance of the Bessel functions is related to the fact that they appear as 
solutions of the field theory equations in various situations. In particular, in spherical and 
cylindrical coordinates the radial parts of the solutions for the scalar, fermionic, and electro- 
magnetic wave equations on background of the Minkowski spacetime are expressed in terms of 
these functions. Another important class of special functions is the so-called Legendre associated 
functions (see, for instance, [21] 122] ). These functions can be considered as generalizations of 
the Bessel functions: in the limit of large values of the degree when the argument is close to 
unity they reduce to the Bessel functions. The associated Legendre functions arise naturally in 
many mathematical and physical applications. In particular, they appear as solutions of physical 
field equations on background of constant curvature spaces (see, for instance, [UEJG3]) and the 
above-mentioned limit corresponds to the limit when the curvature radius of the bulk goes to 
infinity. The eigenfunctions in braneworld models with de Sitter and anti-de Sitter branes are 
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also expressed in terms of the Legendre functions (see [21]). Motivated by this, in the present 
paper, by making use of the generalized Abel-Plana formula, we obtain a summation formula for 
the series over the zeros of the associated Legendre function of the first kind with respect to the 
degree. In particular, this type of series appear in the evaluation of expectation values for phys- 
ical observables bilinear in the operator of a quantum field on background of constant curvature 
spaces in the presence of boundaries. As in the case of the other Abel-Plana-type formulae, 
previously considered in the literature, the formula discussed here presents the sum of the series 
over the zeros of the associated Legendre function in the form of the sum of two integrals. In 
boundary-value problems the first one corresponds to the situation when the boundary is absent 
and the second one presents the part induced by the boundary. For a large class of functions 
the latter is rapidly convergent and, in particular, is useful for the numerical evaluations of the 
corresponding physical characteristics. 

We have organized the paper as follows. In the next section, by specifying the functions 
in the generalized Abel-Plana formula we derive a formula for the summation of series over 
zeros of the associated Legendre function with respect to the degree. In section [3l special 
cases of this summation formula are considered. First, as a partial check we show that as 
a special case the standard Abel-Plana formula is obtained. Then we show that from the 
summation formula discussed in section [21 as a limiting case the formula is obtained for the 
summation of the series over the zeros of the Bessel function, previously derived in [6]. A 
physical application is given in section (4[ where the positive frequency Wightman function 
for a scalar field is evaluated inside a spherical boundary on background of a negative constant 
curvature space. It is assumed that the field obeys Dirichlet boundary condition on the spherical 
shell. The use of the summation formula from section [2] allows us to extract from the vacuum 
expectation value the part corresponding to the geometry without boundaries and to present 
the part induced by the spherical shell in terms of an integral, which is rapidly convergent in 
the coincidence limit for points away from the boundary. The main results of the paper are 
summarized in section In appendix [A] we show that the zeros of the associated Legendre 
function of the first kind with respect to the degree are simple and real, and the asymptotic 
form for large zeros is discussed. In appendix iBl asymptotic formulae for the associated Legendre 
functions are considered for large values of the degree. These formulae are used in section [2] to 
obtain the constraints imposed on the function appearing in the summation formula. 

2 Summation formula 

In this section we derive a summation formula for the series over zeros of the associated Legendre 
function of the first kind, -P^_ 1 / 2 (' u )' w ith respect to the degree, assuming that u > 1 and \i ^ 
(in this paper the definition of the associated Legendre functions follows that given in |22j). For 
given values u and \i this function has an infinity of real zeros. We will denote the positive zeros 
arranged in ascending order of magnitude as z^ : 

1-1/2^) = °' * = 1,2,.... (1) 

These zeros are functions of the parameters u and fi: z^ = Zk(u,fi). Note that one has 

^tz-1/2 ( n ) = ■^-iz-llvft) anc ^' nence > ~ z k are zeros of the function -P^_ 1 / 2 (' u ) as well. In 
appendix [A] we show that the zeros are simple and under the conditions specified above the 
function . ,Au) has no zeros which are not real. 

A summation formula for the series over z\~ can be obtained by making use of the generalized 
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Abel-Plana formula [B] (see also, [3]): 

rb 



f pb 2 ra+ioo 

lim^p.v. / dxf(x) - R[f(z),g(z)] \ = - / dz [g 

MO ° I J a ) 1 Ja-ioD 



(z) + a(z)f(z)}, 



(2) 



where cr(z) = sgn(Imz), the functions f(z) and g(z) are meromorphic for a ^ x ^ 6 in the 
complex plane z = x + and p. v. stands for the principal value of the integral. In formula ([2]) 
we have defined 



R[f(z),g(z)]=7ri 



ERes g(z) + V] Res 

2=z 9,fe z=z f.k 



(3) 



with zj 5 /% and z 9) fc being the positions of the poles of the functions f(z) and g(z) in the strip 
a < x < b. 

The functions f{z) and g(z) in formula ([2]) we choose in the form 



/(*) 
9{z) 



sinh(7r z)h(z) 
h(z) 



-tflTY 



TliP . ,„ fix 

iz — 1/2 V - 



|cos[tt(// + ^)]Qf z _ 1/2 W + cos[vr(/i - iz)]Q^ i2 ,_ 1/2 (V) j 



(4) 



where Qf 2 _i/2( n ) * s ^ ne associated Legendre function of the second kind and the function h(z) is 
meromorphic for a ^ Re z ^ b. By using relation (|70p between the associated Legendre functions 
given in appendix [Bj for the combination appearing on the left hand-side of formula ([2]) one 
finds 

g(z) ± /(z) = — jj cos[vr(^ T w)]Qj te _i /2 (« 



(5) 



With the functions (|4|) the expression for R[f(z),g(z)} takes the form 



R[f(z),g(z)] = 2j2 

k 

with the notation 

r[/i(z)] = 



C-l/ 2 ( n ) 



d z P; 



iz-X/2 



(U) 



cos[7r(/i + iz)]h(z) 



+ 2e- ifin r[h(z)], 



(6) 



2 = Z fc 



(z)iz-X/2 



h(z) 



■ cos[7r(/u — <r(z)iz)]/i(z' 



1 \ - 



Res 



fc.Imz^ fc =0 



*=*fc,fcli^_ 1/2 (tt) 



5^ cos[vr(^ + Uz)]Q? iz _ 1/2 (u) 



l=± 



(7) 



In formula ([7]), z^ ^ are the positions of the poles for the function h(z). 

In terms of the function h(z) the conditions for the generalized Abel-Plana formula ([2]) to 
be valid take the form 



hm / dz—p -^-cos[vr(/i=Fzz)]/i(z) 

'Ja±iw iz-X/2\ U > 
b±ioo Q» Ju) 

dz—^r, cos[7r(/i =F 2Z]Jfi(z) 



lim 



^-1/2 W 



0, 



0. 



(8) 



By using the asymptotic formulae for the associated Legendre functions given in appendix IB"! it 
can be seen that these conditions are satisfied if the function h(z) is restricted to the constraint 



\h(z)\ < e(x)e cm , z = x + iy, \z\ 



oo, 



(9) 
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uniformly in any finite interval of x, where c < 2, e(x)e wx — ► for x — > +00, and 77 is defined 
by the relation 

it = cosh 7/. (10) 

Substituting the functions ([!]) into formula ^ and by taking into account relations ([5]), 
©, we obtain that for a function h(z) meromorphic in the half-plane Re 2; ^ a and satisfying 
condition ([9]), the following formula takes place 



r n 

1 — ™ 



flK Zk,u)h(zk) — p.v. / dx smh(-Kx)h(x) + r[h(z)} 



k=m 



■ pa+ioo (J 1 /r,(u) 



dz — -p — — cos[ir(n - a(z)iz)]h(z), (11) 

Z7T Ja-ioo iz-l/2\ U ) 



where and in what follows the notation 



T ^( z ' u ) = 1/2 \ \ cos Km + **)] ( 12 ) 

° z - r iz-l/2\ U > 



u 



is used. On the left-hand side of formula (jlip . z m —\ < a < z m , z n < b < z n+ \ and in the 
definition of r the summation goes over the poles Zh : k m the strip a < Rez < b. 
Note that from the Wronskian relation for the associated Legendre functions one has 



n » ( v _ e^Tjiz + /z + 1/2) _ 

Vu-ipW- (u 2_i)r(i Z _ M+ i/2)ft,j^_ z ~ Zk - 



Now, by taking into account the formula 

r(iz + ^ + i/2) = ^ |r(iz-^ + i/2)|- 2 

V{iz- n + l/2) cos[vr(/i + iz)\ ' 1 ' 

for the gamma function, the factor T^(zk,u) in (jlip can also be written in the form 

ite^Tiiz - fi+ 1/2) I — 2 



T^(z k ,u) 



Taking the limit a — ► 0, from (jlip one obtains that for a function h(z) meromorphic in the 
half-plane Re z ^ and satisfying the condition © the following formula takes place 

y^r^(^ fc ,u)/t(^ fc ) = —-p.v. / sinh(7rx)/i(x) - r[/i(z)] 
fc=i 2 ^ 

-7T / dx S ' f , co S [ir(n + x)][h(xe™/ 2 ) + h(xe-™/ 2 )}. (15) 

If the function /i(z) has poles on the positive real axis, it is assumed that the first integral 
on the right-hand side converges in the sense of the principal value. From the derivation of 
(|15p it follows that this formula may be extended to the case of some functions h(z) having 
branch-points on the imaginary axis, for example, having the form h(z) = hi(z)/(z 2 + c 2 ) 1//2 , 
where h\(z) is a meromorphic function. This type of function appears in the physical example 
discussed in section [U Special cases of formula (|15p with examples are considered in the next 
section. 



5 



Formula (|15|) can be generalized for a class of functions h{z) having purely imaginary poles 
at the points z = ±iyk, Uk > 0, k = 1, 2, . . ., and at the origin z = yo = 0. Let function h(z) 
satisfy the condition 

h(z) = -h(ze~ m ) + o((z - Ofc)" 1 ), z -> a k , a k = 0, iy fe . (16) 

Now, in the limit a — > the right hand side of ([lip can be presented in the form 

i E (/ E / c K) *«) V , 1 "jt!" ) cos[7r(p " aiz)Mz) > (17) 

plus the sum of the integrals along the straight segments (±z(?/&_i + p),±i(yk — p)) of the 
imaginary axis between the poles. In ([T7|) . C p {ak) denotes the right half of the circle with radius 
p and with the center at the point described in the positive direction. Similarly, 7+ and 7" 
are upper and lower halves of the semicircle in the right half-plane with radius p and with the 
center at the point z = 0, described in the positive direction with respect to this point. In the 
limit p — > the sum of the integrals along the straight segments of the imaginary axis gives the 
principal value of the last integral on the right-hand side of (j!5p . Further, in the terms of (|17|) 
with a = — we introduce a new integration variable z' = ze m . By using the relation (I16p the 
expression (|17p is presented in the form 

^ f Qti _i io( u ) 1 

- (l-W2)Res p" < \ ^s[n(p-iz)]h(z)\ (18) 

* k =o,iy k Z ~ ak 1 VVJ W J 

plus the part which vanishes in the limit p — > 0. As a result, formula (|15p is extended for 
functions having purely imaginary poles and satisfying condition (|16p . For this, on the right- 
hand side of (Q2D we have to add the sum of residues (fTHj) at these poles and take the principal 
value of the second integral on the right-hand side. The latter exists due to condition (|16p . Note 
that for functions having the form h(z) = F(z)P^ z _ l ^ 2 (u) the left-hand side of (I15p is zero and 
from this formula we obtain a formula relating the integrals involving the Legendre associated 
functions. 



3 Special cases 

Here we will consider special cases of the summation formula f 1 1 5 [) . First let us consider the case 
p = —1/2. The corresponding associated Legendre functions have the form 

In this case one has zt = Trk/rj. Introducing a new function F(x) in accordance with the relation 
F^tjx/tt) = sinh(Trx)h(x), and assuming that this function is analytic in the right half-plane, 
from formula (1151) we find the Abel-Plana formula in the standard form: 



J2m = -lm+ rdxF(x)+t r d X F{ix \ nx F{ . ix) . (20) 

k=i 2 Jo Jo e -1 

Note that the first term on the right-hand side of this formula comes from the residue term with 
<J k = in (US]). 
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In the case \i = 1/2 for the corresponding associated Legendre functions we have the expres 
sions 



' ^-i/ 2 ( coshr ?) = * 



o (21) 

7r y / smhr/ ' " Sz ~ 1 / 2V 1 y 2 -y/sinhr? 

The zeros z^ now have the form Zk = Tt(k + 1/2) /r, and for functions F(z) analytic in the right 
half-plane from formula (|15p we obtain the Abel-Plana formula in the form useful for fermionic 
field calculations (see, for instance, [HE]): 



£>(fc + l/2) 



k=l 



dx F{x) 



00 , F(ix) - F(-ix) 
dx- 



-,2-kx 



+ 1 



(22) 



As a next special case let us consider the formula for the summation over the zeros of the 
function P~^_^ / 2 (cosh(r//s)) in the limit when s — * oo. By taking into account the relation (see 
appendix iBj) 



lim 



(23) 



with Ju(ji) being the Bessel function of the first kind, in this limit from (|15p we obtain the 
summation formula for the series over zeros r] = juk, k = 1,2,..., of the Bessel function. In 
order to take this limit we also will need the formulae ()77[) from appendix [B] and the formulae 

m\ 



lim v ii P u M [cosh(x/z/)] 
lim i/'QjJ'^cosh^/i/)] 



e'^K^x), 



(24) 



with Iu(x), Kn(x) being the modified Bessel functions. First we rewrite formula (115j) making 
the replacements z — ► sz, x — > sx, fi — * —fi, in both sides of this formula including the terms in 
r[/i(z)], and we take u = cosh(^/s). In order to take the limit s — > oo for the second integral on 
the right-hand side of the resulting formula, we note that, as it follows from the derivation of 
(fT5|) . the integrand of this integral (with the replacements described above) should be understood 
as the limit 



cos[tt (s x — fx)] ^2 h(sxe lm / 2 ) = lim^ cos[-7r(sx — \i — lise)]h( 



sxe 



(25) 



*=+,- 



Taking the limit s — > oo with the help of formulae (|23p . (|24j l.( )77p . we find the following summation 
formula over the zeros of the Bessel function 



p.v. 

1 

7T 



dxf(x) - n[f(z)] 
Kjx 



dx 



e^7(xe™ /2 ) + e'^fixe-™/ 2 ) 



(26) 



where f(z) = lim^oo e sz ^ ri h(sz/rj), and 
rj[f(z)} = 



i Res 

^Im z h<k >0 

-ir y Res 



H?\z) 



V Res 



—Res 

2 z=o 



(27) 
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This formula is a special case of the result derived in [5] (see also, [3]). 

Now let us consider two important special cases of (|15p corresponding to [i = —I, h(z) = 
H(z)/ cosh(7rz) and fj, = —1 — 1/2, h(z) = H(z)/ sinh(7rz) with I = 0, 1,2, . . .. The associated 
Legendre functions with these values of the order appear as radial solutions of the equations for 
various fields on background of constant curvature spaces in cylindrical and spherical coordinates. 
Let the function H(z) is meromorphic in the half-plane Rez ^ and satisfy the condition 



\H(z)\ < £H(x)e cvy , z = x + iy, \z\ 



CO, 



(28) 



uniformly in any finite interval of x > 0, where c < 2, £h( x ) —> for x — > +oo. Then from the 
results of section [2] it follows that the formula 



E 

k=l 



z=z k 



1 

5 P.v. 



dx tanh 1_5 (vrx)i?(x) - r s [H(z) 



— / dx 

2ir 



Q 



-1-8/2 
x-1/2 



(«) 



H{xe™l 2 ) + (-l) d H(xe-^'%(29) 



■ni/2\ 



x-1/2 



takes place, where 5 = 0, 1. In this formula we have introduced the notation 



E 



fc,Im Zh,k¥=0 

+ E 

k,Imz h fc=0 



Res 

Z=Zh,k 



o*(z) 



Q 



-1-8/2 
-a{z)iz-l/2 



U 



r-z-^M 



-H(z) 



Res 

z=z h , k 



H(z] 



(30) 



Adding to the right-hand side of formula (|29p the term 



{-1) & Y, (l-W2)Res 



Q 



-1-8/2 
-iz-1/2 



(«) 



P 



-1-5/2 



iz-1/2 



(u) 



-H(z) 



(31) 



with H(z) obeying the condition H(z) = — (—l) s H(ze~ m ) + o((z — ak)~ 1 ) for z — > a^, we obtain 
the extension of this formula to the case when the function H(z) has poles at the points 0, ±yk- 
Prom (|29p . as an example when the series is summarized in closed form one has 



E 



-.-1-1/2, \ 

Qiz-x/2 W z 2n cos(az) 



id z P+)' 2 \u) i^ + c^ 



\m+n 



z=z k 



ml 



d m 
dx r 



7T 



d_ 

2 m+2 y c Q c 
-1-1/2 



, c 2n-l e - aC) 



Qx-i/2 ( u ) x 2n cosh(ax) 



(x + c] 



m+1 



(32) 



where a < 2rj, c > 0, with m ^ and ^ n ^ m being integers. The last term on the right-hand 
side of this formula comes from the residue at the pole = ic. As a next example, we take in 
formula (1291) with 5 = 1 the function 



H(z) 



J2n—v 



J v {az) 



( z 2 +c 2)a/2 ' 



(33) 
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where a, b, c are positive constants and n is a non-negative integer. This function is analytic in 
the right half-plane and satisfies the condition ()28[) if a + b < 2rj, 2n < a + v. By taking into 
account that (|33|) is an even function of z, from (|29p we find 



(^ + c 2 )«/ 2 5,^-^(11)1^ 2i (z 2 + c 2 )*/ 2 

4 Wightman function inside a spherical boundary in a constant 
curvature space 

In this section we consider a physical application of the summation formula derived in section 
[21 Namely, we will evaluate the positive frequency Wightman function for a scalar field and 
the vacuum expectation value of the field squared inside a spherical shell in a constant negative 
curvature space assuming that the field obeys the Dirichlet boundary condition on the shell (for 
quantum effects on background of constant curvature spaces see, for instance, [H [23] and 
references therein). 

Consider a quantum scalar field <p(x) with the curvature coupling parameter £ on background 
of the space with constant negative curvature described by the line element 

ds 2 = dt 2 - a 2 [dr 2 + sinh 2 r{d8 2 + sin 2 Ode/) 2 )] , (35) 

where a is a constant which is related to the non-zero components of the Ricci tensor and the 
Ricci scalar by the relations 

12 ^ ^ 
R x = R 2 = R z = R = (36) 

a z a z 

The field equation has the form 

'ViV 1 + M 2 + £R) <p(x) = 0, (37) 



where M is the mass of the field quanta. 

We are interested in quantum effects induced by the presence of a spherical shell with radius 
r = ro, on which the field obeys Dirichlet boundary condition: (p(x)\ r=ro = 0. This boundary 
condition modifies the spectrum of the zero-point fluctuations compared with the case of free 
space and changes the physical properties of the vacuum. Among the most important char- 
acteristics of the vacuum are the expectation values of quantities bilinear in the field operator 
such as the field squared and the energy-momentum tensor. These expectation values are ob- 
tained from two-point functions in the coincidence limit. As a two-point function here we will 
consider the positive frequency Wightman function W(x,x') = {0\(p(x)<p(x')\0) , where |0) is the 
amplitude for the vacuum state. This function also determines the response of Unruh-De Witt 
type particle detectors [23]. Expanding the field operator over the complete set {(p a (x), (x)} 
of classical solutions to the field equation satisfying the boundary condition, the Wightman 
function is presented as the mode-sum 

W(x,x') = J2<P*(x)<PUx'), (38) 
a 

where a is a set of quantum numbers specifying the solution. 
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By the symmetry of the problem under consideration, the eigenfunctions for the scalar field 
can be presented in the form 

< Pa (x) = Z{r)Y lm (9,<l>)e- iut , (39) 

where Yi m (6,(j)) are standard spherical harmonics, I = 0,1,2,..., —l^m^l. From the field 
equation (|37l) we obtain the equation for the radial function Z(r): 



1 d ( 2 dZ \ 

smh r— — + 



r 2 2 \ 2 1(1 + 1) 
[uj -m eS )a - , 

sinn r 



Z = 0, (40) 



sinh 2 rdr\ dr 

where we have introduced the effective mass defined by 

m e 2 ff = M 2 — 6£/a 2 . (41) 

In the region inside the spherical shell the solution of equation (|40|) . finite at r = 0, is expressed 
in terms of the associated Legendre function of the first kind and the eigenfunctions have the 
form 

~P ^ f cosh. 7**") 

<p a (x) = C a lZ - l/ L— Y lm (9, fie-™, (42) 

V smh r 

with the notation 

z 2 = (uj 2 - m 2 cS )a 2 - 1. (43) 

From the boundary condition on the spherical shell we find that the eigenvalues for z are 
solutions of the equation 

p-i- 1 / / 2 (coshro)=0, (44) 

and, hence, z = Zk, k = 1, 2, . . ., in the notations of section[2j The corresponding eigenfrequencies 
are found to be 

oj 2 k = oj 2 {z k ) = (z 2 + 1 - 60/a 2 + M 2 . (45) 

Hence, the set a of the quantum numbers is specified to a = (I, m, k). 

The coefficient C a in (|42p is determined from the orthonormalization condition 

d 3 x y/\g\y a (x)w* a ,{x) = (46) 

where the integration goes over the region inside the spherical shell. Substituting the eigen- 
functions (|42jl into (I46p . by taking into account the integration formula (j65H and the boundary 
condition, one finds 

c - 2 = a 3^1 (n 2 _ l)^p-i-^(no)^p-i- 1 / 2 2 (^)U = , ft , u=uo , (47) 

where and in the discussion below we use the notations 

u = coshr, uq = coshro. (48) 
By using the Wronskian relation (|13p . the formula for the normalization coefficient is written as 

, fc r(^ + / + l)Qr;-_Y 2 2 (n )e^/ 2 )- 

Cl = _ ; _ 1/2 • (49) 

a 3 u(z k )T(iz k - l)d z P iz _ x j 2 (u )\ z=Zk 

Note that the ratio of the gamma functions in this formula can also be presented in the form 

T(iz k + 1 + 1) 2 cos[ir(iz - I - 1/2)] 
T{ lZk - I) = lF{lZk + 1 + 1)1 n ' (50) 
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Substituting the eigenfunctions into the mode-sum formula ([55]) and using the addition 
theorem for the spherical harmonics, for the Wightman function one finds 

= ' f; < 2 ' + p^e'^t^T^ + I + 1)|' 

47r z a d f— ' Vsmhrsmhr' f— f 

(=0 fc=l 

— i— 1/2 -2-1/2 e -Mz k )At 

x T-*-i/ 2 (z fc , «o)P, fc _ 1/2 (cosh r )P. 2fc _ 1/2 (cosh r ) ^ , (51) 

where At = t — t' and T fM (z,u) is defined by relation (fl2j) . In ([51~]) . p(cos7) is the Legendre 
polynomial and 

cos 7 = cos 6 cos 0' + sin 6 sin #' cos(</> — 4>'). (52) 

As the expressions for the zeros z^ are not explicitly known, formula (|5ip for the Wightman 
function is not convenient. In addition, the terms in the sum are highly oscillatory for large 
values of quantum numbers. 

For the further evaluation of the Wightman function we apply to the series over k the 
summation formula (|15p taking in this formula 

-iuj(z) At 

w-1/2 V""" ' > x w-1/2 V— * ' I u ( z j 

The corresponding conditions are met if r + r' + At /a < 2rg. In particular, this is the case in 
the coincidence limit t = t' for the region under consideration. For the function (|53p the part of 
the integral on the right-hand side of formula (|15|) over the region (0, xm) vanishes and for the 
Wightman function one finds 



h(z) = z\T(iz + I + 1)\ 2 P-J^i 2 (cosh r)F;i- 1 / 2 2 (cosh r'f—— -. (53) 



47r ° ^ Vsmhrsmhr' A u r(x — f) 



^x-i/f^o) , 1/a , 1/2 cosh(Jx 2 - x 2 M At/a) 
x *_( - P^-// 2 (cosh t)P-^ (cosh /) V , (54) 

Vs-1/2 ^ V X X M 



where we have defined 

x M = VM 2 a 2 + 1 - 6f. (55) 



In formula (|54p . the first term on the right-hand side is given by 

W (x,x>) = 1 y(!+^W e -» rdxxsinh(vrx) 
tlo vsinhrsinhr' Jo 

p—ico(x)At 

x\T(ix + I + l)| 2 P^(coshr)P;^ 1 7^coshr0^ r . (56) 

This function does not depend on the sphere radius and is the Wightman function for a scalar 
field in background spacetime described by the line element (|35p when boundaries are absent. 
This can also be seen by the direct evaluation. Indeed, when boundaries are absent the eigen- 
functions are still given by formula (|42p . where now the spectrum for z is continuous. In this 
case the corresponding part on the right of the orthonormalization condition (|46p should be 
understood as the Dirac delta function. In the case z = z' the normalization integral diverges 
and, hence, the main contribution comes from large values r. By using the asymptotic formulae 
for the associated Legendre functions for large values of the argument, we can see that 

2 



T(iz) 



T{l + l + iz) 



5(z - z'). (57) 
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By using this result for the normalization coefficient in the case when boundaries are absent one 
finds 



1 



V2: 



tua° 



T{l + l + iz) 



Viz) 



and the eigenfunctions have the form (see also, [H [26] ) 



<p a (x) 



T(l + l + iz) 



T(iz) 



P. 



-1-1/2 
:-l/2 



(cosh r) 



V 2ua 3 sinh r 
(EH 



-Y, 



lm \ 



~iu)t 



(58) 



(59) 



for the corresponding Wightman 



Substituting these eigenfunctions into the mode-sum 
function we find the formula which coincides with (|56p . 

The case of a spherical boundary in the Minkowski spacetime is obtained in the limit a — > oo, 
with fixed or = R. In this limit one has xm = aM. Introducing a new integration variable 
y = x/a, using the formulae (|24p and the asymptotic formula for the gamma function for large 
values of the argument, we find 



W^ M \x,: 



w, 



(M) 



(x,X f ) - ^ 



(2Z + l)P,(cos 7 ) 



1=0 



dyy 



M 



'ii+i/2(Ry)ii+i/2(R'y) 



Att 2 ^/RR' 
Ki+i/ 2 (Roy) cosh(Vy 2 - M 2 At) 



(60) 



This formula gives the the positive frequency Wightman function inside a spherical shell with 
radius i?o i n the Minkowski bulk and is a special case of the general formula given in the first 
paper of [14J for a scalar field with Robin boundary conditions in arbitrary number of spatial 
dimensions. 

Having the Wightman function (|54p. we can evaluate the vacuum expectation value of the 
field squared taking the coincidence limit of the argument. Of course, this limit is divergent 
and some renormalization procedure is necessary. Here the important point is that for points 
outside the spherical shell the local geometry is the same as for the case of without boundaries 
and, hence, the structure of the divergences is the same as well. This is also directly seen 
from formula (154j) . where the second term on the right-hand side is finite in the coincidence 
limit. Since in formula (|54p we have already explicitly subtracted the boundary-free part, the 
renormalization is reduced to that for the geometry without boundaries. In this way for the 
renormalized vacuum expectation value of the field squared one has 



'0,ren 



J(1+1/2)tt 



{21 + 1) 



Att 2 c 



sinhr 



dx x 



T(x + l + l) 



Qji-J/^coshro) fciJ^CcoBhr) 



T(x-l) P - 



1-1/2 
-1/2 



(cosh ro) 



(61) 



'M 



where the first term on the right-hand side is the corresponding quantity in the constant nega- 
tive curvature space without boundaries and the second one is induced by the presence of the 
spherical shell. For large values x, the integrand in (|6ip behaves as e~^ r °~ r ^ x /(2x sinhr) and the 
integral is exponentially convergent at the upper limit for strictly interior points. 

For r — ► one has Pji^ 2 (cosh r) ~ (r/2)' +1 / 2 /r(/ + 3/2), and in the boundary induced 
part at the sphere center the I = term contributes only: 



'O.ren 



1 

2^c 



dx ■ 



x [x 



)_ 

?2xro _ i 



-M 



-1/2 



0. 



(62) 
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where we have used formulae (|19p . Note that for a conformally coupled field the boundary 
induced part in (|62p coincides with the corresponding quantity for the sphere with radius aro 
in the Minkowski bulk. 

5 Conclusion 

The associated Legendre functions are an important class of special functions that appear in a 
wide range of problems of mathematical physics. In the present paper, specifying the functions 
in the generalized Abel-Plana formula in the form (J3]), we have derived summation formula 
(|15p for the series over the zeros of the associated Legendre function -P^_ 1 / 2 (' u ) wr th respect to 
the degree. This formula is valid for functions h(z) meromorphic in the right half-plane and 
obeying condition ([9]). Using formula (|15p . the difference between the sum over the zeros of the 
associated Legendre function and the corresponding integral is presented in terms of an integral 
involving the Legendre associated functions with real values of the degree plus residue terms. 
For a large class of functions h(z) this integral converges exponentially fast and, in particular, 
is useful for numerical calculations. Frequently used two standard forms of the Abel-Plana 
formula are obtained as special cases of formula (1151) with [i = — 1/2 and [i = 1/2 and for an 
analytic function h(z). Applying the summation formula for the series over the zeros of the 
function Pj 1 z _ 1 ^ 2 (ucosh(r]/s)) and taking the limit s — > oo we have obtained formula (I26p for 
the summation of the series over zeros of the Bessel function. The latter is a special case of the 
formula, previously derived in [6]. Further, we specify the summation formula for two special 
cases of the order [i = — I and [i = —1 — 1/2 with / being a non-negative integer and give examples 
of the application of this formula. The associated Legendre functions with these values of the 
order arise as solutions of the wave equation on background of constant curvature spaces in 
cylindrical and spherical coordinates. 

In section H] we consider a physical application of the summation formula. Namely, for a 
quantum scalar field we evaluate the positive frequency Wightman function and the vacuum ex- 
pectation value of the field squared inside a spherical shell in a constant negative curvature space 
assuming that the field obeys the Dirichlet boundary condition on the shell. In spherical coordi- 

^ 1/2 

nates the radial part of the corresponding eigenfunctions contains the function P iz _-±_ % (coshr) 
and the eigenfrequencies are expressed in terms of the zeros by relation (|45p . As a result, 
the mode-sum for the Wightman function includes the summation over these zeros. For the 
evaluation of the corresponding series we apply summation formula (|15p with the function h(z) 
given by (|53[) . The term with the first integral on the right-hand side of formula (|15p corresponds 
to the Wightman function for the constant curvature space without boundaries and the term 
with the second integral is induced by the spherical boundary. For points away from the shell 
the latter is finite in the coincidence limit and can be directly used for the evaluation of the 
boundary induced part in the vacuum expectation value of the field squared. The latter is given 
by the second term on the right-hand side of formula (|6ip . The renormalization is necessary 
for the boundary-free part only and this procedure is the same as that in quantum field theory 
without boundaries. 

On the physical example considered we have demonstrated the advantages for the application 
of the Abel-Plana-type formulae in the evaluation of the expectation values of local physical 
observables in the presence of boundaries. For the summation of the corresponding mode-sums 
the explicit form of the eigenmodes is not necessary and the part corresponding to the boundary- 
free space is explicitly extracted. Further, the boundary induced part is presented in the form of 
an integral which rapidly converges and is finite in the coincidence limit for points away from the 
boundary. In this way the renormalization procedure for local physical observables is reduced to 
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that in quantum field theory without boundaries. Note that methods for the evaluation of global 
characteristics of the vacuum, such as total Casimir energy, in problems where the eigenmodes 
are given implicitly as zeros of a given function, are described in references |27| . 
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A On zeros of the function , ,Ju) 

iz— 1/2 v / 

In this appendix we show that the zeros z = Zk are simple and real. By making use of the 
differential equation for the associated Legendre functions it can be seen that the following 
integration formula takes place 

Taking the limit v' — > v and applying the Lopital's rule for the right-hand side, from this formula 
we find 

/*. = u - .^^Mi^M-gwww + const (64) 

J 1v + 1 

By taking into account the relation P f ^ iz _ 1 ^ 2 (u) = ^z-1/2^)' we see that ^ or rea ^ z one ^ as 
[Pi z _ 1/2 (u)} 2 = \P£_ 1/2 (tt)\ 3 . Hence, from formula §4$ we find 

£ dv\P? z _ 1/2 (v)\ 2 = {[d z Pt z _ 1/2 (u)]d u Pt z _ 1/2 (u) - P»_ 1/2 (u)d z d u P? z _ 1/2 (u)} . (65) 

Here we have taken into account that for u — > 1 one has Pi z _ij 2 { u ) ~ ( M — 1)~ M and, hence, 
lim^^i Pi z _i/ 2 ( u ) = for fj, < 0. From formula ([65]) it follows that [d z Pj 1 z _ 1 ^ 2 (u)] z=Zk / 0, and, 
hence, the zeros Zk are simple. 

Now let us show that under the conditions u > 1 and /j, ^ all zeros of the function 
^2-1/2 ( M ) are rea ^' Suppose that z = A is a zero of P^ z _i/ 2 (u) which is not real. As the function 
P z _i/ 2 {u) has no real zeros (see, for instance, |28|). A is not a pure imaginary. If A* is the 
complex conjugate to A, then it is also a zero of Pi z _i/ 2 { u )i because Plx*_i/ 2 (v) = [^A-i/2^)]*' 
As a result, from formula f|63|) we find 

f*C-i/ 2 (^A-i/2(«) = 0- (66) 

We have obtained a contradiction, since the integrand on the left hand-side is positive. Hence 
the number A cannot exist and the function P. M , ,„(u) has no zeros which are not real. 

tz— 1/2 V ' 

From the asymptotic formula (|7ip for the function Pj* z _i/ 2 (u) (see appendix IBl below) we 
obtain the asymptotic expression for large zeros: 

z k ~ (irk - nfi/2 - tt/4) /r/. (67) 

Note that this result can also be obtained by taking into account that for large values z from 
(|23|) one has Pr^, 2 (cosh (77)) ~ z~^ J^(r]z) and using the asymptotic form for the zeros of the 
Bessel function (see, for instance, [22]). 
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B Asymptotics of the associated Legendre functions 

In this appendix we consider asymptotic expressions for the associated Legendre functions for 
large values of the degree. As a starting point we use the formula 

rf 1 11 4- 7 4- i A (1 — p^ 2r i)^ 

<#_ 1/3 (cosM) = ypae^ r(i + .) e (^/2), F W 2 + ^ V2 + * + m; i + *; e" 2?? ). (68) 

Using the linear transformation formula 15.3.4 from [22] for the hyper geometric function, the 
expression for the function (cosh 7/) is presented in the form 

g^ 1/2 (cosh^) = 0^ V2 sinh/ (1/2 + M ' 1/2 " M; 1 + Z] 1/(1 " e2r?)) - (69) 

Now, by using the result that for large \c\ one has F(a,b;c;z) = 1 + 0(l/|c|), from (|6T?j) the 
asymptotic formula for the function (5^_ 1//2 (coshry) is obtained for large values \z\. The corre- 
sponding formula for the function P^ 1 _ 1 , 2 (coshry) is obtained by using the relation 

7re lM7r sm(jiz)P^ ^^(coshr?) = cos[7r(z — n)]Q^ z ^(coshr?) — cos[7r(z + fJ-)]Q^ _ 1 ^ 2 (coshr/). (70) 
In this way we obtain the following formulae 

PI yi*- 1 / 2 

P M _ 1/2 (coshr/) ~ y--^===sin(r7y-i7?x + 7r/i/2 + 7r/4), 
fir 

^-i/2( coshr ?) ~ y 2 v^ee^ exp [ -7?x ~ ~ ~ ^Z 4 ^' ( 71 ) 

in the limit y — ► +oo, z = x + %, and the formulae 

T M-l/2 

e i/2 (oosh,) ~ yfe-^le— ™, (72) 

in the limit cc — > +oo. 

Now let us consider the asymptotics of the functions P i ~^ 1 ^ 2 (cosh(ry/z^)), Q^ i/ _ 1 , 2 (cosh(r]/i')) 
as v — * +oo. These asymptotics are obtained in the way similar to that used in [25] for formulae 
(|24p . Our starting point is the formula 

F^ 1/2 (cosh(ry/2,)) = ^^F(l/2 - iu, 1/2 + iv; 1 + M ; - sinh 2 (r//2z,)), (73) 

relating the associated Legendre function to the hypergeometric function. From the definition 
of the hypergeometric function it is not difficult to see that 

lim F(l/2-zV,l/2 + w/;l + /x;-sinh 2 (77/2i/)) =T(1 + (74) 

v -^+00 

Combining (|73h and (174)) we obtain formula (|23f) . The corresponding formula for the functions 
<5^f i/ _ 1/ / 2 ( cosn ( r ?/ zy )) are obtained by making use of the relation 

-sin( / uvr)e^Q7' 1 u Ju) = T ^ ±lv ~ ^ + l / 2 \ P ^ ( u ) - P7\ ._(«), (75) 
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and formula (|23p . In this way we find 

or in the equivalent form 

Jim^i^e^Q-^^CcoBhCT//!/)) = -^e^H^irj), 
^u^QZl^coBh^/v)) = ^-e^H^ir,), (77) 

(1 2) 

where ' (rj) are the Hankel functions. 
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